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A  NEW  COMPUTATIONAL  CAPABILITY  FOR  RAMJET  PROJECTILES 


Sukumar  R.  Chakravarthy  and  Uriel  C.  Goldberg 

Rockwell  International  Science  Center 
Thousand  Oaks,  California 

Abstract 

/ 

A  state-of-the-art  computer  code  has  been  applied  to  calculate  internal  and  external 

flows  for  a  ramjet  projectile  under  turbulent,  supersonic  conditions.  A  multizone  compu- 
WJ 

tational  approach ^ has  beei*  used,  and  its  advantages  are  demonstrated.  Two  turbulence 
modeling  approaches  were  tried:  1)  The  Baldwin-Lomax  model  applied  throughout  the  flow 
field;  2)  The  Baldwin-Lomax  model  applied  outside  of  backflow  regions  and  a  dedicated 
algebraic  Jfc-e  backflow  model  applied  within  these  regions.  Comparison  with  experimental 
data  shows  a  significant  improvement  in  the  code’s  predictional  capability  when  the  latter 
approach  is  used.  (  ) 

1.  INTRODUCTION 

/ 

Solid  fuel  ramjet  technology  is  currently  of  high  interest.  Initial  tests  for  spin- 
stabilized  projectiles  demonstrated^fche  feasibility,  of  the  technology.  However,  many  funda¬ 
mental  issues  must  still  be  resolved  m~ordef  to  properly  exploit  this  emerging  technology. 
For  example,  combustion  efficiency  as  observed  in  flight  is  much  lower  than  preliminarily 
predicted.  This  may  be  due  to  the  spin  of  the  projectile.  On  the  other  hand,  the  pre¬ 
liminary  predictions  may  have  been  inaccurate  because  the  internal  flow  and  combustion 
processes  of  hollow  ramjet  projectiles  are  not  completely  understood.  Without  a  thorough 
knowledge  of  such  processes,  scaling  laws  cannot  be  established  to  design  new  projectiles . 
that  may  be  fin  or  spin-stabilized.  A  series  of  “cold  flow”  wind  tunnel  tests  (no  combus¬ 
tion)  have  been  conducted,  and  earlier  numerical  computations  have  not  shown  consistent 
agreement  with  the  wind  tunnel  measurements.  Wind  tunnel  tests  are  very  expensive 
and  new  ramjet  models  can  not  be  established  solely  on  the  basis  of  empirical  data  bases 
from  previous  flight  histories  and  wind  tunnel  tests.  New  accurate  and  reliable  numerical 
methods  must  be  established  to  support  the  ramjet  program. 

Many  numerical  methods  have  artificial  viscosity  (sometimes  called  numerical  smooth¬ 
ing  or  dissipation)  controlled  by  selecting  various  coefficients.  The  presence  of  these  terms 
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and  their  wide  range  of  possible  parameters  are  highly  undesirable.  Often  individual  so¬ 
lutions  must  be  carefully  fine  tuned  to  achieve  acceptable/stable  solutions.  These  dissipa¬ 
tive  terms  become  very  troublesome  in  high  speed  flow  problems  when  turbulence  models 
(which  also  have  several  adjustable  constants)  must  be  applied.  The  inherent  complex 
geometries  associated  with  solid  fuel  ramjets,  together  with  these  other  problems,  present 
a  very  difficult  simulation  problem.  The  number  of  free  parameters  within  the  numerical 
solution  must  be  reduced  to  a  minimum. 

It  is  clear  that  numerical  methods  with  more  exact  modeling  (fewer  adjustable  or  em¬ 
pirical  terms)  must  be  developed  for  the  ramjet  problem.  Recently  a  new  class  of  numeri¬ 
cal  methods  has  been  developed  that  are  more  exact  and  reliable.  These  approaches  have 
been  termed  total  variation  diminishing  (TVD)  formulations.  Such  numerical  methods 
are  under  investigation  by  a  few  researchers  for  cases  of  inviscid  flow  and  two-dimensional 
or  simple  geometries  (ramjet  problems  need  three  dimensions,  have  nozzles,  inlets,  and 
other  irregular  geometrical  features  such  as  flame  holders).  The  present  authors  have 
published  these  types  of  solutions,  and  have  since  produced  a  series  of  codes  where  vis¬ 
cous/turbulent/separated  flow  processes  within  complex  geometries  are  treated.  As  will  be 
shown  in  this  paper,  the  presence  or  absence  of  an  adequate  turbulence  model  dedicated 
to  the  treatment  of  separated  flow  regions  makes  all  the  difference  between  satisfactory 
and  poor  predictional  capability  of  the  state-of-the-art  computer  code  used  in  the  present 
work. 


2.  COMPUTATIONAL  APPROACH 

The  computational  fluid  dynamics  (CFD)  approach  is  to  use  state-of-the-art  nu¬ 
merical  algorithms  based  on  the  TVD  formulation  along  with  an  implicit,  factored  time¬ 
stepping  algorithm  in  a  zonal  grid  framework.  The  TVD  formulation  contributes  to  a 
very  reliable  (no  numerical  oscillations,  no  “numerical”  user  knobs  such  as  dissipation  co¬ 
efficients)  and  accurate  methodology.  The  solution  procedure  using  implicit  schemes  also 
contributes  to  computational  efficiency.  The  use  of  zonal  grids  simplifies  the  construction 
of  computational  grids,  allows  complex  topologies  to  be  treated  routinely  and  permits  more 
options  in  working  with  limited  computer  main  core  memory,  etc. 

The  TVD  formulation  is  used  to  discretize  the  hyperbolic  part  of  the  Navier-Stokes 
equations.  The  diffusion  components  of  these  equations  are  treated  with  central  difference 
approximations  for  second  derivative  terms  in  each  direction.  The  cross-derivative  terms 
are  discretized  in  a  novel  fashion  to  augment  the  diagonal  dominance  of  the  resulting  set 
of  difference  equations. 


2.1  Computational  Grid  —  A  Zoned  Approach 

The  problem  of  computing  internal/external  ramjet  projectile  flow  fields  is  compli¬ 
cated  by  the  complex  geometry  and  the  complex  physics  involved.  Resolving  the  physics 
in  complex  geometries  could  necessitate  the  use  of  dense  grids  for  which  the  required  data 
storage  might  easily  outstrip  available  in-core  computer  memory.  In  the  present  approach, 
a  zonal  grid  is  used  to  treat  the  ramjet  intemal/extemal  flow  geometries.  This  allows  con¬ 
siderable  flexibility  in  responding  to  the  complexities  of  geometry,  physics,  and  computer 
resources. 

In  the  zonal  approach,  the  computational  method  and  the  computer  program  are 
constructed  in  such  a  manner  that  each  zone  may  be  considered  as  a  more-or-less  inde¬ 
pendent  module  with  the  zones  interacting  with  each  other  after  or  before  the  information 
corresponding  to  each  zone  is  updated  one  cycle. 

Complex  geometries  are  gridded  more  easily  by  zonal  approaches.  The  complete 
physical  domain  to  be  considered  can  be  broken  up  into  simpler  topological  subdivisions. 
Even  a  single  topological  subunit  can  be  divided  into  many  zones.  Selective  and  adaptive 
grid  refinement  can  be  employed  in. some  of  these  zones  to  increase  accuracy  and  resolution. 

Zonal  approaches  are  also  one  way  of  reducing  the  required  in-core  computer  memory. 
In  such  approaches,  it  would  be  sufficient  to  have  enough  main  core  storage  to  fit  the  needs 
of  a  single  zone. 

The  usefulness  of  a  particular  grid  can  be  verified  by  repeating  the  calculations  on 
twice  finer  or  coarser  a  grid  and  estimating  the  truncation  error  from  the  two  sets  of 
data.  The  zonal  approach  permits  such  an  estimation  to  take  place  in  a  restricted  region 
of  the  entire  domain  with  an  attendant  saving  in  computational  resources  (both  time 
and  memory).  Zonal  approaches  also  permit  the  easy  implementation  of  local  refinement 
techniques.  Specific  regions  of  the  flow  field  may  be  selected  for  additional  refinement  and 
can  be  treated  as  a  new  zone. 

The  multigrid  strategy  will  be  employed  at  a  later  stage  to  enhance  computational 
efficiency  even  further.  The  multigrid  method  involves  the  use  of  a  hierarchy  of  crude  and 
fine  grids  in  a  cyclical  pattern,  the  aim  being  to  reduce  the  residue  on  the  finest  grid  by 
performing  calculations  on  coarser  grids.  The  computer  program  organization  required 
for  multigrid  techniques  is  identical  to  that  required  for  zonal  approaches  and,  thereby, 
the  zonal  approach  enables  a  natural  transition  from  zonal  methods  to  zonal  multigrid 
methods. 


In  each  zone,  any  approach  to  generating  the  computational  grid  may  be  used:  al¬ 
gebraic  methods,  differential  equations  methods,  etc.  In  the  latter  category,  elliptic  or 
hyperbolic  methods  may  be  used,  as  applicable.  The  multisurface  technique  is  representa¬ 
tive  of  the  former  type  of  grid  generation  method. 

This  subsection  on  grid  generation  is  now  concluded  by  actually  illustrating  the  mul¬ 
tizone  grid  treatment  for  a  ramjet  projectile  geometry.  Figure  1  shows  the  division  of  the 
complete  geometry  into  6  zones.  Figure  2  shows  the  composite  6-zone  grid  used  for  the 
calculations  and  consisting  of  the  following  meshes  in  Zones  1  through  6:  30  x  35,  7  x  17, 
50  x  30,  30  x  17,  15  x  25,  and  45  x  20.  The  grids  have  been  generated  in  such  a  manner  that 
the  horizontal  families  of  mesh  lines  within  the  projectile  are  continuous  between  zones 
and  the  vertical  families  are  continuous  in  the  zones  external  to  the  ramjet  projectile.  This 
type  of  grid  is  natural  for  the  problem  at  hand  but  the  computational  methodology  permits 
a  wider  class  of  patched  zonal  grids:  neither  family  of  grid  lines  need  be  continuous  across 


2.2  Computational  Algorithm 

The  computational  discretization  technique  must  be  reliable  and  robust  if  it  is  to 
successfully  capture  the  complex  physics  of  ramjet  internal  and  external  flow  fields.  The 
TVD  formulation  for  convection  terms  (the  hyperbolic  part  of  the  time-dependent  Navier- 
Stokes  equations),  along  with  a  careful  treatment  of  diffusion  terms,  offers  itself  as  the  best 
alternative. 

TVD  formulations  are  relatively  new  but  well  proven  and  documented1-4.  By  design, 
they  have  some  desirable  features:  they  are  highly  accurate  (up  to  third-order  accuracy 
for  steady-state  solutions  to  hyperbolic  equations);  they  avoid  numerical  oscillations  which 
could  arise  in  high  gradient  regions,  near  shock  waves,  etc.;  they  avoid  the  formation  of 
unphysical  expansion  shock  waves;  they  contribute  to  the  diagonal  dominance  of  the  finite 
difference  equations.  They  are  completely  defined;  there  are  no  “numerical”  dissipation 
and  other  coefficients  to  be  tuned  case  by  case.  They  have  successfully  been  applied  to 
compute  a  very  wide  range  of  inviscid  and  viscous  flows  without  recourse  to  such  “tuning 
knobs”. 

Proper  treatment  of  diffusion  terms  that  arise  in  the  Navier-Stokes  equations,  in  the 
equations  of  turbulence  modeling,  species  equations  of  combustion  chemistry,  etc.,  are  also 
important  in  the  construction  of  reliable  numerical  methods.  Unidirectional  second  deriva¬ 
tive  terms  are  conveniently  handled  correctly  by  using  central  difference  approximations. 


Cross  derivatives  are  conventionally  treated  in  cavalier  fashion.  In  the  present  work,  a 
careful  cross  derivative  treatment  which  augments  diagonal  dominance  is  used2.  Once 
again,  this  is  a  novel  but  proven  technique  which  will  contribute  much  to  the  reliability 
and  robustness  of  the  numerical  procedure. 

Conventional  approximate  factorization  methods  are  used  to  handle  multidimensional 
implicit  formulations  and  this  is  the  preferred  approach  in  the  present  development.  How¬ 
ever,  the  diagonal  dominance  of  TVD  schemes  (with  properly  discretized  cross  derivative 
terms)  makes  it  possible  to  use  relaxation  methods5. 

Many  techniques  can  be  applied  to  implement  an  efficient  code  based  on  the  above 
methods.  The  equations  of  turbulence  modeling  will  be  decoupled  from  the  Reynolds- 
averaged  Navier-Stokes  equations  in  the  implicit  inversion  procedure.  This  decoupling 
will  have  the  effect  of  reducing  the  magnitude  of  the  computational  linear  algebra  problem 
of  inverting  a  block  tridiagonal  system  of  equations  and  will  not  have  a  significant  effect 
on  the  time  accuracy  of  the  overall  numerical  procedure.  The  entire  code  is  “vectorized” 
(which  involves  a  restructuring  of  FORTRAN  loops,  eliminating  recursive  and  conditioned 
statements  in  those  loops,  etc.)  to  a  very  high  degree  and  will  be,  therefore,  very  suitable 
for  high  speed  vector  computers  such  as  the  CRAY-XMP.  Other  techniques  which  may  be 
used  selectively  include:  the  use  of  spatially  varying  time  steps  to  accelerate  the  solution 
to  time-asymptotic  steady  states  when  possible;  the  use  of  “diagonalization”  approaches 
in  which  the  left-hand  side  of  the  implicit  approximately  factored  form  of  the  difference 
equations  is  preconditioned  and  the  resulting  equations  are  solved  as  a  set  of  scalar  equa¬ 
tions.  (This  approach  requires  the  viscous  terms  to  be  placed  only  on  the  right-hand  side 
of  the  equations.  It  also  has  a  side  effect  on  the  conservation  property  but  only  in  the 
transient  stages  and  not  when  steady  state  is  reached.) 


auations  of  Motion 


The  compressible,  Reynolds-averaged  Navier-Stokes  equations  for  2D/axisymmetric 
flow  are  written  in  the  following  conservation  form,  where  the  dependent  variables  u,  v,  and 
e  are  mass-averaged,  with  e  being  the  specific  total  internal  energy,  T  the  temperature,  p 
and  p  being  mean  density  and  pressure,  respectively,  and  t  being  time: 


dW  dF_  dG  ( G  _  H\ 
dt  +  dx  +  dy  +  \  y  y  ) 
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In  Eq.  (1),  the  laminar  and  eddy  viscosities,  p  and  pt,  are  implicitly  divided  by  the 
reference  Reynolds  number.  The  equations  used  for  the  Euler  calculations  are  obtained 
from  Eq.  (1)  by  setting  both  laminar  and  eddy  viscosities  to  zero. 


In  all  the  calculations,  the  air  was  assumed  to  be  a  perfect  gas,  satisfying  the  equation 
of  state 

P  =  pRT.  (2) 


The  following  power  law  was  used  to  relate  molecular  viscosity  to  temperature6: 

p/po  =  (T/T0)n  (3) 


where  p0  =  0.1716  mP,  T0  =  491.6°R,  and  n  =  0.64874. 
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The  laminar  and  turbulent  Prandtl  numbers,  Pr  and  Prt,  were  assumed  constant 
with  values  of  0.72  and  0.9,  respectively.  The  ratio  of  specific  heats,  7,  was  also  assumed 
constant  and  equal  to  1.4.  Cv  and  Cp  are  specific  heat  capacities  at  constant  volume  and 
constant  pressure,  respectively. 

For  time-invariant  grids,  in  the  £-7  computational  plane,  Eq.  (la)  is  transformed  into 
the  finite  volume  conservation  law  form  represented  by2 


I 

& 


/- 

.V 


dW  1 
dr  +  Area 


(yvF-xnG)t 


,  r,  ™  G  H 
+  (-y^F  +  X{G)V  - 

y  yj 


=  0 


(4) 


where  £  and  77  are  the  new  independent  variables  and  x^,  xv,  y$,  and  yn  are  the  four 
transformation  coefficients  obtained  numerically  from  the  mapping  procedure.  The  “Area” 
in  Eq.  (4)  denotes  the  area  of  the  finite  volume  cell  under  consideration  at  the  time  of 
discretization  of  the  equations,  r  denotes  transformed  time. 


3.  TURBULENCE  MODELING 

Ramjet  internal  flows  involve  large  regions  of  recirculatory  flow,  induced  both  by  shock 
waves  and  by  sharp  geometrical  discontinuities.  Most  existing  turbulence  models  either 
do  not  treat  such  regions  or  do  so  in  an  ad  hoc  fashion  which  is  frequently  inadequate.  A 
notably  different  approach  is  the  use  of  a  full  Reynolds  stress  closure  model,  involving  the 
solution  of  five  coupled  partial  differential  equations  (for  two-dimensional  flows)  for  the 
three  normal  stresses,  the  shear  stress,  and  the  length  scale.  These  must  be  supplemented 
by  a  wall  function  to  provide  turbulence  quantities  across  viscous  regions  adjacent  to  solid 
surfaces.  Such  a  wall  function  is  usually  some  form  of  the  law -of- the- wall,  which,  according 
to  experimental  observations,  does  not  apply  to  separated  flows.  Thus,  an  expensive  and 
time  consuming  computation  of  the  Reynolds  stresses  is  coupled  with  a  questionable  near- 
wall  formulation. 

In  an  attempt  to  improve  predictional  capability  of  separated  flows  using  current 
Reynolds-averaged  Navier-Stokes  codes,  a  new  and  simple  turbulence  model  has  recently 
been  introduced7,  some  details  of  which  given  below. 

3.1  Highlights  of  the  Turbulence  Model 

The  new  turbulence  model  is  based  on  experimental  observations  of  separated  tur¬ 
bulent  flows.  The  model  prescribes  turbulence  kinetic  energy  (k)  and  its  dissipation  (e) 
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analytically  inside  separation  bubbles.  A  Gaussian  variation  of  k  normal  to  walls  is  as¬ 
sumed.  The  length  scale  of  turbulence  within  bubbles  is  proportional  to  the  local  distance 
from  the  wall  to  the  edge  of  the  viscous  sublayer,  which  is  located  outside  the  backflow 
region,  as  shown  in  Fig.  3.  The  latter  feature  is  a  basic  assumption  of  the  model. 

The  velocity  scale  is  the  local  maximum  Reynolds  stress,  which  typically  occurs  around 
the  middle  of  the  boundary  layer,  well  outside  the  bubble.  This  scale  must  be  supplied  by 
a  turbulence  model  which  is  used  beyond  separated  regions. 

The  main  equations  of  the  backflow  model  are  given  in  the  Appendix.  A  simple 
formula  for  eddy  viscosity  distribution  within  separation  bubbles  results,  and  is  used  to 
supply  eddy  viscosity  for  the  Reynolds-averaged  equations  when  performing  the  calcula¬ 
tions  inside  the  bubbles.  Outside,  another  turbulence  model  (e.g.,  Baldwin-Lomax  or  k-e) 
supplies  the  values  of  eddy  viscosity.  In  all  the  calculations  presented  here,  the  Baldwin- 
Lomax  model  was  applied  outside  backflow  regions.  In  order  to  mimic  the  influence  of 
the  large  eddies  residing  outside  separation  zones,  the  following  averaging  procedure  is 
performed  when  switching  from  the  separation  model  to  the  one  used  beyond  backflows: 
outside  backflows,  normal  to  walls,  the  eddy  viscosity  is  taken  to  be  the  average  of  the  one 
supplied  by  the  model  used  beyond  backflows  and  that  given  by  the  separation  model  with 
/  =  A  +  B  (see  Appendix);  downstream  of  reattachment,  the  history  effect  of  the  large 
eddies  is  taken  into  account  by  averaging  the  locally  predicted  eddy  viscosity  with  the  one 
provided  by  the  separation  model  at  the  location  within  backflows  where  the  skin  friction 
attains  its  maximum  negative  value. 

The  mutual  influence  of  multiple  wails  on  the  eddy  viscosity  is  taken  into  account  by 
using  an  inverse  averaging  procedure: 


(£."=,  v»;f 


where  N  is  the  number  of  walls,  rii  is  the  local  normal  distance  from  the  ith  wall,  and  fiti 
is  the  eddy  viscosity  due  to  the  presence  of  the  »th  wall  by  itself. 

Further  details  about  the  new  model  and  its  performance  are  given  in  Refs.  7-9. 
It  should  be  noted  that  the  model  is  not  of  an  eddy  viscosity  type;  it  can  be  used  in 
conjunction  with  Reynolds  stress  models. 
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The  computer  code  described  previously  has  been  applied  to  calculate  both  the  in¬ 
ternal  and  the  external  flow  fields  for  the  ramjet  configuration  shown  in  Fig.  1  under  the 
following  conditions:  Afoo  =  4.03,  Reoo  =  7  x  107  /meter.  All  walls  were  assumed  to  be 
adiabatic. 

The  code  was  applied  in  three  modes:  1.  inviscid  (Euler);  2.  turbulent,  using  the 
Baldwin-Lomax  (B-L)  model  everywhere;  3.  turbulent,  using  the  B-L  model  outside  back- 
flow  regions  and  the  new  backflow  model  within  these  regions. 

Figure  4  compares  internal  wall  pressure  predictions  under  these  three  modes  with 
experimental  data  from  BRL10  for  a  1.7  in.  injector  diameter.  Both  the  inviscid  prediction 
and  the  one  using  the  B-L  model  by  itself  are  poor,  grossly  underpredicting  the  pressure 
rise  along  the  cavity  wall.  The  prediction  which  incorporates  the  backflow  model,  however, 
is  much  improved  throughout  the  length  of  the  internal  wall. 

Figure  5  shows  skin  friction  distribution  along  the  internal  wall  as  resulting  from  the 
Mode  3  run.  The  extent  of  the  reversed  flow  regions  is  clearly  seen.  Since  the  Mode  3 
option  proved  superior  to  the  others,  all  subsequent  calculations  were  performed  using  this 
mode. 

Figure  6  is  a  pressure  contour  plot  for  the  ramjet  projectile  with  1.534  in.  nozzle 
diameter  and  1.9  in.  injector  diameter.  The  details  of  the  flow  field  are  captured  very 
sharply,  a  fact  which  is  fully  attributable  to  the  TVD  scheme.  Note  the  oblique  shock  in 
the  inlet  section.  In  Fig.  7,  the  inner  wall  pressure  profile  prediction  is  compared  with  the 
experimental  data10.  Overall,  the  agreement  is  very  good,  although  the  pressure  upstream 
of  the  injector  is  somewhat  overpredicted.  This  may  be  the  result  of  an  incorrectly  pre¬ 
dicted  shape  of  the  leading  edge  of  the  separation  bubble  (see  next  figure)  which  would 
produce  an  oblique  shock  of  larger  strength  than  the  observed  one.  Figure  8  shows  the 
corresponding  skin  friction  distribution,  indicating  that  most  of  the  flow  upstream  of  the 
injector  as  well  as  the  entire  cavity  flow  are  detached  from  the  wall. 

The  next  calculation  was  done  for  a  geometry  consisting  of  a  1.1  in.  nozzle  diameter 
with  a  1.7  in.  injector  diameter.  The  pressure  contour  plot  of  Fig.  9  clearly  indicates  a 
normal  shock  just  upstream  of  th'  cowl  inlet,  a  result  of  internal  flow  restriction  due  to 
the  reduced  nozzle  throat  size.  This  renders  the  flow  subsonic  all  the  way  to  downstream 
of  the  nozzle  throat.  Figure  10  shows  inner  wall  pressure  distribution  for  this  case.  Two 
data  sets  are  indicated:  one  corresponding  to  the  1.1  in.  nozzle  diameter  used  in  the 
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calculations,  the  other  to  a  1.2  in.  nozzle  diameter  case.  The  predicted  pressure  agrees 
with  the  latter  rather  than  with  the  former.  Assuming  that  the  data  are  reliable,  this 
would  indicate  a  nozzle  wall  displacement  thickness  distribution  which  is  predicted  to  be 
thinner  than  the  actual  one.  This  could  be  the  result  of  inadequate  grid  resolution  on  the 
nozzle  wall,  deficiency  of  the  turbulence  model,  or  both.  Further  studies  will  be  performed 
to  clarify  this  point. 

Figure  11  shows  the  sharpness  at  which  the  normal  shock  is  captured,  a  direct  conse¬ 
quence  of  the  TVD  scheme.  The  figure  also  shows  the  outer  wall  pressure  distribution  in 
the  vicinity  of  the  cowl,  where  the  initial  rise  due  to  the  oblique  shock  and  the  subsequent 
drop  due  to  the  wall  curvature  are  clearly  seen.  In  Fig.  12,  the  inner  wall  skin  friction 
distribution  is  shown.  This  time  only  the  cavity  flow  is  detached  from  the  wall. 

Finally,  it  was  found  important  to  indicate  the  sensitivity  of  the  measured  data  to 
small  changes  in  nozzle  throat  diameter.  Figure  13  shows  two  sets  of  data  for  the  1.7  in. 
injector  diameter  case:  one  with  a  1.534  in.  nozzle  throat  diameter,  the  other  with  a 
1.564  in.  diameter.  This  1.9%  change  produces  two  completely  different  internal  flow 
fields,  as  seen  by  the  pressure  distributions.  Also  included  in  the  figure  is  a  calculation 
corresponding  to  the  1.534  in.  nozzle.  The  pressure  distribution  lies  between  the  two  data 
sets  but  closer  to  the  data  corresponding  to  the  1.564  in.  nozzle  case.  This  sensitive  range 
of  nozzle  diameters  will  be  studied  further  and  reported  in  future  work. 

5.  CONCLUSIONS 

A  state-of-the-art  computer  code  has  been  applied  to  calculate  internal  and  external 
flows  for  a  ramjet  projectile  under  turbulent,  supersonic  conditions,  using  a  multizone 
technique.  Two  types  of  turbulence  modeling  approaches  were  used:  1.  the  Baldwin- 
Lomax  model  by  itself;  2.  the  Baldwin-Lomax  model  applied  outside  of  backflow  regions 
in  conjunction  with  a  dedicated  backflow  model  applied  within  such  regions.  Comparison 
with  experimental  data  showed  a  significant  improvement  in  predictional  capability  when 
the  backflow  model  was  invoked.  Since  the  latter  model  is  a  simple  algebraic  one,  it 
should  be  found  attractive  for  applications  which  involve  recirculating  flow  regions.  Further 
work  is  planned  to  assess  the  influence  of  grid  resolution  and  turbulence  modeling  on  the 
predictional  capability  of  the  code.  The  sensitivity  of  the  internal  flow  to  small  changes  in 
nozzle  diameter  will  also  be  studied. 
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APPENDIX 


A  SUMMARY  OF  BACKFLOW 
TURBULENCE  MODEL  EQUATIONS 

Following  are  the  most  important  equations  on  which  the  separation  turbulence  model 
is  based. 

k  and  e  within  the  separation  bubble  are  given  by 

pk/Pikh  =  [e+/(e+  -  1)]  [l  -  «-«•/"*>’] 

=  G(*,n)  ,  0  <  n  <  nk, 

(p/Ph)3,2e  =  [KG(s,nj\3/2 /nv  ,  0  <  n  <  nk, 

where 

<t>  =  0.5,  0  =  kv/kh  =  1  +  [(n,/n»)J  -  l]^/(e*  -  1), 
pvkv/pw  =  u]/yfci,nl  =  20Cj/4  +  nj  , 
n*  =  nutfuw  ,  Cp  =  0.09  ,  C *  =  0.7, 

and  «,  =  (-«'«' )m«  =  Here,  is  the  maximum  Rey¬ 

nolds  stress,  assumed  to  correspond  to  the  maximum  normal-to-wall  mean  velocity  gradi¬ 
ent.  vt,m  is  the  value  of  eddy  viscosity  where  the  aforementioned  maximum  gradient  occurs 
(see  Fig.  3),  and  it  must  be  supplied  by  the  turbulence  model  which  is  applied  external  to 
the  separation  bubble. 

The  length  scale  of  turbulence  within  the  bubble  is  given  by  L  =  n„//?3/2. 

Finally,  eddy  viscosity  is  given  by 

=  yj^  [f(n/nb)/( Sv^/?2)]  n*vGl/2(s,  n)  , 

0  <  n  <  n4 


where  f(n/nh)  =  A(n/nb)  +  B,A  =  -(C;/2)9/5  ,  B  =  (C;/2)3/5  -  A. 
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This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  consents/answers  to  the  iteas/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Nuaber _ Date  of  Report _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Consent  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ _ _ 
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- - - - -  t 

S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
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etc?  If  so,  please  elaborate. _ _ 


6.  General  Consents.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 
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